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The amount of water vapor in the atmosphere helps to determine the likelihood that 
severe storms may develop. The concentration of water vapor, though, is highly variable 
in space and time. And yet small changes in water vapor concentration over a short 
period of time or over a short spatial distance can dete rmin e whether a storm may or may 
not develop. Therefore, in order to improve the ability to forecast severe weather such as 
thunderstorms it is important to measure water vapor in the atmosphere with high spatial 
and temporal resolution. 

One of the most attractive research tools for measuring water vapor in the atmosphere 
with high spatial and temporal resolution is a Ra ma n lidar. A Raman lidar consists of a 
laser transmitter, a telescope receiver and optics and electronics for processing optical 
and electronic signals. A laser pulse is emitted into the atmosphere and it interacts with 
molecules in the atmosphere causing them to become excited and to emit, through the 
Raman process, photons of different wavelength than emitted by the laser. The molecule 
that emitted these photons can be identified based on the wavelength of the photons 
emitted. This is the way that a Raman lidar identifies water vapor molecules in the 
atmosphere. 

One of the great challenges in Raman lidar measurements has been to make useful 
daytime measurements of the water vapor profile under bright daytime conditions. In this 
first of two papers, we describe the instrumentation and analysis of the first documented 
Raman lidar that is able to measure water vapor in the daytime with sufficient quality to 
permit the study of developing storm systems. 
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Abstract 

The NASA/GSFC Scanning Raman Lidar (SRL) participated in the International H 2 0 Project (IHOP) 
that occurred in May and June, 2002 in the Midwestern part of the U. S. The SRL received extensive 
optical modifications prior to and during the IHOP campaign that added new measurement capabilities and 
enabled unprecedented daytime water vapor measurements by a Raman Lidar system. Improvements were 
also realized in nighttime upper tropospheric water vapor measurements. The other new measurements that 
were added to the SRL for the IHOP deployment included rotational Raman temperature, depolarization, 
cloud liquid water and cirrus cloud ice water content. In this first of two parts, the details of the operational 
configuration of the SRL during IHOP are provided along with a description of the analysis and calibration 
procedures for water vapor mixing ratio, aerosol depolarization and cirrus cloud extinction to backscatter 
ratio. For the first time, a Raman water vapor lidar calibration is performed taking full account of the 
temperature sensitivity of Raman water vapor and nitrogen scattering. Part two presents case studies that 
permit the daytime and nighttime error statistics to be quantified. 
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2 Introduction 


Hie International H 2 O Project was the largest meteorological field campaign ever held in the United States 
[Wechworth et. al., 2004]. It occurred in the mid-western United States between May 13 - June 25, 2002. 
The goal of IHOP was to measure convective storm systems with sufficient detail to permit quantitative 
precipitation and convection initiation forecasting to be improved. During IHOP, numerous intensive op- 
erations periods were declared focussing on boundary layer evolution, drylines, bores, nocturnal jets, and 
other mesoscale events. IHOP included participants from numerous U.S. and foreign government agencies 
as well as universities. 

The instru me ntation used during IHOP included seven research aircraft carrying three water vapor li- 
dars and one wind lidar, mobile radar systems for storm-chasing, and a ground-based site in the western 
panhandle of Oklahoma that included the NASA/Goddard Space Flight Center (GSFC) Scanning Raman 
Lidar (SRL), Goddard Laboratory for Winds (GLOW) molecular wind lidar [Gentry et. al., 2004], and the 
Holographic Airborne Rotating Lidar Instrument Experiment (HARLIE) [Schwemmer et. al., 2004] scan- 
ning aerosol lidar. Other instruments that were located at the western Oklahoma site that came to be 
known as "Homestead" were Vaisala RS-80 G1 5H and SnowWhite radiosonde systems [Wang et. al., 2003], 
cloud and aerosol radar, sodar. Atmospheric Emitted Radiance Interferometer (AERI) [Feltz et. al., 1998], 
SuomiNet Global Positioning System (GPS) [Ware et. al., 2000] and an Aerosol Robotic Network (AERONET) 
[Holben et. al., 1998] sun photometer. 

This first of two papers provides a detailed description of the experimental configuration of the NASA/GSFC 
SRL during its participation in IHOP. This configuration included measurements of water vapor mixing ra- 
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tio, aerosol backscatter, extinction and depolarization, rotational Raman temperature measurements per- 
formed in the UV [Di Girolamo et al., 2004], cloud liquid water [Russo et. al, 2004], and cirrus cloud 
particle size and ice water content [Wang et al., 2004]. All of these measurements were made using a 
single output wavelength of 354.7 nm demonstrating the feasibility of providing these measurements in an 
automated and eye-safe system [Di Girolamo et. al., 2004], Particular attention will be given here to the 
instrumentation and analysis techniques that led to improved daytime water vapor measurement capability. 
With these improvements, it is now possible to use a single experimental configuration for quantifying wa- 
ter vapor variation in the boundary layer during the daytime and extending into the upper troposphere at 
night. The paper is structured as follows. A brief history of the use of the narrow-band, narrow field-of- 
view technique for Raman water vapor lidar measurements is provided as context for the SRL water vapor 
measurements in IHOP. The techniques for processing the raw data and of calculating water vapor mix- 
ing ratio, aerosol scattering ratio, aerosol depolarization and cirrus cloud extinction to backscatter ratio are 
then presented. In part two of this paper, comparisons of the SRL water vapor measurements with those of 
other instruments at IHOP will be presented along with daytime and nighttime case studies that demonstrate 
measurement capability and permit the error statistics to be quantified. 

3 The Pursuit of Non-Solar Blind Daytime Raman Water Vapor Lidar 
Measurements 

During IHOP, the SRL demonstrated greatly improved non-solar-blind daytime Raman water vapor lidar 
measurements over any previously demonstrated. The combination of technologies and techniques t ha t per- 
mitted this were 1) a large-pulse tripled Nd:YAG laser, 2) narrow field of view telescope, 3) narrow spectral 
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band detection, 4) fast photomultiplier detectors and 5) a combination of analog and photon counting elec- 
tronics, which permit the full intensity Raman signals to be sampled. While none of these elements is new, 
the combination had not been fully exploited previously for Raman water vapor lidar. Furthermore, the 

I; . . . • 

path toward this solution was not direct and occurred over more than 20 years of research activity at various 
research centers. 

3.1 The first efforts in Raman water vapor lidar at NASA/GSFC 

The first Raman water vapor lidar measurements were made nearly simultaneously in 1969 in separate re- 

i 

search efforts by S. H. Melfi [Melfi et. al., 1969] and J. Cooney [Cooney, 1970]. Development of Raman 

lidar technology for measurement of water vapor was started at NASA/GSFC in 1980 when Dr. Melfi 

j 

joined the center. The first Raman water vapor lidar system developed at GSFC was based on a Quan- 
taRay DCR1-A tripled Nd:YAG laser (354.7 nm) operating at 10 Hz and providing 150 mj per pulse in 

| 

the UV. The laser had a single, flash-lamp pumped, amplifier and used harmonic crystals for converting 
the 1.064 pm fundamental wavelength to the doubled and tripled wavelengths of 532 nm and 354.7 nm, 

| 

respectively. The crystal housings were not temperature stabilized necessitating tilt-tuning of the crystals 
to maximize the output energy every two minutes during extended data acquisition periods. The telescope 
I used was based on a 1.5 m searchlight mirror that was built in Japan and used on a German surface ship in 

World War II. The mirror was found in a military equipment storage facility at Chesapeake Beach, Virginia, 
moved to GSFC, re-coated and placed into service for use in Raman lidar. A Cassagrain configuration 
telescope was constructed using this mirror so that the detector package could be installed under the very 
large telescope. The telescope had relatively poor optical quality so that the received signals were measured 
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using a field of view of 2-5 milliradians. Furthermore, due to uncertainty about whether the substrate ma- 
terial had been annealed, no hole was cut in the mirror so that the return signal were transmitted through 
the uncoated substrate. A simple detector package consisting of beamsplitters, interference filters and two 
photomultiplier tubes measuring the Raman signals from water vapor (407.5 run) and nitrogen (386.7 nm) 
was used. The interference filters were manufactured by Barr Associates and possessed approximately 
15-20% transmission with bandwidths of 5-10 nm. The data system used LeCroy 8837F 8-bit, 32-MHz 
analog-to-digital (AD) transient recorders and was controlled by a Digital Equipment Computer PDP-11 
running the RT-1 1 operating system. This system was known as the Marine Boundary Layer Raman Water 
Vapor Lidar and was successfully used to chart the evolution of water vapor in the boundary layer in 1985 
[Melfi and Whiteman, 1985]. The system was housed in a mobile trailer and was first used in a field de- 
ployment for the Cooperative Huntsville Meteorological Experiment (COHMEX) in Huntsville, Alabama 
in 1986 [Dodge et. al., 1986] [Fulton and Heymsfield, 1991] [Straka and Anderson, 1993]. 

The nighttime limitation of this early Raman water vapor lidar system was clear but there was no real 
progress toward the development of daytime Raman water vapor lidar measurements at GSFC throughout 
the 1980s despite the feet that the solar blind technique of making daytime Raman water vapor measure- 
ments was being demonstrated successfully by others [Renault and Capitini, 1988]. During the 1980s, up- 
grades to the Marine Boundary Layer Raman Water Vapor Lidar included replacing the searchlight-based 
telescope with an astronomical quality, 0.75 m aperture system and implementing dual photomultiplier de- 
tection channels in a low altitude/high altitude configuration for water vapor, nitrogen and Rayleigh-Mie 
measurements. The new telescope, although possessing just 1/4 of the collecting area of the former one, 
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significantly improved the signal-to-noise ratio of the measurements due to the greatly reduced losses (both 
scattering and transmission) and improved focussing quality. The data acquisition system was also up- 
graded to add photon counting (PC) electronics (Joerger S3 200 MHz counters) to the pre-existing analog 
detectors. The first published measurements of water vapor and aerosol structure of cold and warm frontal 
passages were made using this system [Melfi et. al., 1989], although the measurements were still limited to 
nighttime. 


By the late 1980s, most of the necessary technology was in place to attempt a non-solar blind measure- 
ment of daytime water vapor except that interference filter technology still did not permit the fabrication 
of the high transmission, narrow bandwidth optics that are now available. Furthermore, frequent ground 
loop and impedance matc hing problems associated with the simultaneous use of AD and PC measurements 
encouraged testing the idea of using photon counting electronics exclusively as a replacement for the com- 
bined AD and PC system that was in use previously. Two photomultipliers tubes, both operating in the 
photon counting mode, were used for each measured wavelength. The near range channel was attenuated 
significantly to permit it to be photon-counted while the far-range channel received the full intensity sig- 
nal. This concept produced excellent measurements during nighttime tests performed using the Marine 
Boundary Layer Water Vapor Udar in 1987. Due to funding limitations in the late 1980s, development 
work on Raman lidar at GSFC ceased for a period of time and this mobile Raman water vapor lidar was 
converted into the first mobile stratospheric ozone lidar: STROZ-LITE (Stratospheric Ozone Lidar Trailer 
Experiment) [McGee et. al., 1991]. 
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3.2 Early Non-Solar Blind Daytime Raman Water Vapor Lidar Activities 


The first published measurements of daytime water vapor mixing ratio using non-solar blind Raman lidar 
were performed in 1991 in northern Germany [Ansmann et. al., 1992] using the narrow field-of-view, nar- 
row spectral band technique being discussed here. The measurements were made using a high power XeCl 
laser (>67.5W), 300 MHz photon counting electronics, 0.8 m telescope operating at 0.1 - 0.4 milliradian 
field-of-view, and grating polychromator providing approximately 0.3 nm bandpass. An additional edge 
filter was used in the water vapor channel to increase the suppression of the Rayleigh-Mie signal. Due to 
the low optical efficiency of the system, no attenuation of the signal was required during daytime measure- 
ments to permit the signal to be recorded using only photon counting electronics (private communication, 
Drs. Albert Ans mann and Ulla Wandinger). Measurements with 10-50% random error were made with 
this system to an altitude of ~1 .7 km using 15 minutes averaging, and ~200 m spatial smoothing. 

The need within the atmospheric radiation community to improve the measurements of water vapor and 
other parameters to support atmospheric radiative transfer modeling led to technology development within 
the U.S. Department of Energy (DOE) of water vapor measurement systems with improved accuracy and 
precision [Goldsmith et. al., 1998], The DOE Instrument Development Program (IDP) began in 1989 and 
provided partial fiinding for the development of the first implementation of the NASA/GSFC Scanning 
Raman Lidar (SRL). This work was carried out jointly between NASA/GSFC and DOE Sandia National 
Laboratory (SNL) in a research effort that resulted in the construction of two Raman lidar systems: one 
at GSFC and the other at SNL. The systems shared the same telescopes and used excimer lasers, but key 
aspects of the designs were different to explore multiple techniques for making daytime Raman water vapor 
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lidar measurements. The NASA/GSFC effort pursued the solar-blind approach for daytime water vapor 
measurements and therefore used a fluorine-optimized excimer laser that could operate using either a XeF 
(35 1 nm for nighttime measurements) or a KrF (248 nm for daytime measurements) gas mixture. The SNL 
effort used the narrow-band, narrow field of view approach [Bisson et. al., 1999] and was based on a XeCl 
laser providing 40-50W of output power. This parallel development effort resulted in the first side-by-side 
demonstration of Raman water vapor lidar measurements in 1994 [Goldsmith et.al., 1994], The SNL system 
later demonstrated the concept of using a large aperture telescope at narrow field of view for improved 
daytime measurements [Bisson et. al., 1999] and helped to demonstrate some of the concepts that were 
designed into the CART Raman Lidar, the first automated Raman lidar system [Goldsmith et al., 1998]. 

The success of this technology demonstration effort led to the development of a new Raman lidar at 
the DOE Southern Great Plains (SGP) site in northern Oklahoma [Goldsmith et. al., 1998]. That system, 
referred to now as CARL (Climate Research Facility Raman Lidar), became operational in 1996 and uses 
a NdrYAG laser operating at 9-12W and the narrow-band, narrow-field-of-view approach. CARL’s original 
configuration used only photon counting data acquisition thus requiring aq factor of 10 attenuation of the 
Raman signals during the daytime to prevent photon saturation. Nonetheless, with a 10-minute average and 
75-meter range resolution, water vapor mixing ratio measurements with less than 20% random error to an 
altitude of approximately 3 km were performed. Several years of such measurements are now available, pro- 
viding continuous measurements of water vapor throughout the diurnal cycle [Turner and Goldsmith, 1999], 
In June of 2004, the data acquisition system of CARL was upgraded to include both analog and photon 
counting electronics essentially identical to those in use in the SRL for IHOP permitting the full strength 
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Raman signals to be used during both daytime and nighttime [?]. 

4 The Scanning Raman Lidar (SRL) 

Development of the NASA/GSFC Scanning Raman Lidar began in 1989 under the support of both NASA 
and DOE’s IDP Program. The SRL was first deployed in the field in November 1991 for the Spectral 
Radiance Experiment [Ellingson and Wiscombe, 1996 ] that was sponsored by NASA and DOE and took 
place in Coffeyville, Kansas. At that time the system was based on a Lambda Physik LPX 240iCC excimer 
laser that produced up to 30 watts at 351 nm using a XeF gas mixture and up to 100 Watts at 248 ran using 
a KrF gas mixture. Nighttime measurements were performed using the 351 nm output of the laser and 
daytime measurements [Whiteman et. al., 1994] were made using the 248 nm output and the solar blind 
technique [Renault and Capitini, 1988] whereby solar background is essentially eliminated by absorption 
due to stratospheric ozone. The use of the solar blind technique requires a knowledge of the tropospheric 
ozone profile so that the differential extinction of the H 2 O and N 2 signals by ozone can be accounted for in 
the water mixing ratio calculation. This approach to Raman water vapor measurements permitted photon 
counting data acquisition (100 MHz in the initial implementation of the SRL) to be used exclusively. The 
solar blind measurements of water vapor using the SRL were hampered by the large attenuation of the 
outgoing laser radiation due to tropospheric ozone which caused signal-induced-noise (SIN) problems in 
the PMT detectors in use at the time. The PMTs were upgraded in 1993 to ones with significantly lower 
SIN problems, but still overlap problems prevented an accurate quantification of the ozone profile in the 
lowest 1 km [Whiteman et. al., 1994]. Also, the absorption cross section of ozone does not differ greatly 
between the Raman-shifted returns for N 2 and O 2 when excited by 248 nm [Goldsmith and Fenare, 1994] 
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[Whiteman et. al., 1994]. This degrades the sensitivity of the technique for deriving ozone. Theoretical 
modeling [Goldsmith and Ferrare, 1994] indicates that a longer excitation wavelength of approximately 
26Qnm would be greatly preferred for solar blind measurements. Nonetheless, daytime solar blind water 
vapor mixing ratio measurements were made with the SRL in 1993 to an altitude of 2.5 km using a 20 
min ute average [Whiteman et al., 1994]. 

The narrow spectral band, narrow field-of-view technique for making Raman water vapor lidar measure- 
ments was also explored with the SRL in 1992 using the XeF (351 ran) output of the excimer laser. These 
measurements were performed by expanding the laser beam by a factor of 5 and using a 0.5 milhradian field 
of view. However, the 100-MFIz photon counting data acquisition system, fully sufficient for the original 
solar-blind conception of the SRL, was incapable of handling the high countrates present in non-solar blind 
daytime Raman lidar measurements. Thus, the input signals required attenuation by more than 2 orders 
of magnitude to permit them to be photon-counted. Nonetheless, daytime measurements were performed 
through the boundary layer with 1 hour averaging [Whiteman et. al., 1994], In the early 1990’s therefore we 
had determined that addition of analog electronics would be necessary to optimize daytime measurements 
using the narrow-band, narrow-field-of-view technique. 

It was also clear that, considering a fixed amount of output power, a smaller number of larger laser 
pulses offered significant advantages for daytime operations [Bisson et. al., 1999] over a larger number of 
smaller pulses. Therefore, in 1995, a Nd:YAG laser was added to the SRL in preparation for the TARFOX 
(Tropospheric Aerosol Radiation Field Observational Experiment) field campaign [Ferrare et. al., 2000]. 
Narrow spectral band, narrow field-of-view measurements of water vapor were performed using the SRL 
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during TARJFOX and can be seen at our website (http://ramanlidar.gsfc.nasa.gov), however the 100-MHz 
photon counting data system was still in use requiring that the input signals be attenuated by at least an 
order of magnitude in order to limi t the photon arrival rate. In addition, the interference filters available at 

the time were limited to approximately 25% transmission for a 0.5 nm passband. 

4.1 Additional technology development and the SRL configuration for IHOP 

Drawing on the previous experience gained in Raman water vapor lidar development, in 1998 we began 
construction of the Raman Airborne Spectroscopic Lidar (RASL) under the NASA Instrument Incubator 
Program (IIP). The goal of this instrument was to bring the range of Raman lidar water vapor and aerosol 
measurements to an airborne platform and to take advantage of performance enhancements that are realized 
by measuring downward in the atmosphere versus measuring upward [Whiteman et. al., 2001b]. This new 
development effort included fee purchase of a new data acquisition system from Licel of Berlin, Germany 
that incorporates both analog and photon counting electronics and that circumvented one of fee limitations 
of fee earlier daytime Raman water vapor lidar measurements. The RASL data acquisition system was used 
in fee SRL for fee IHOP field deployment. By the time of IHOP, significantly improved interference filters 
were also available for Raman water vapor measurements. These two enhancements were critical to the 
successful water vapor measurements made by fee SRL during IHOP. 

Many additional modifications were made to the SRL immediately prior to and during fee IHOP field 
campaign. Substantially new wavelength separation optics were implemented on fee large 0.76m F/5.2 
Dall-Kirkham telescope. This provided simultaneous measurements of water vapor mixing ratio, aerosol 
backscatter and extinction, cloud liquid [Whiteman and Melfi, 1999] [Russo et. a], 2004] and ice water 
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[Wang et. al., 2004], and rotational Raman temperature [Di Girolamo et. al., 2004], In addition, a 0.25m 
F/2.5 Newtonian telescope was added in the field during the IHOP experiment for the low altitude Raman 
measurements, accomplished via fiber optic coupling, and for aerosol depolarization measurements using 
free-space coupled optics. The fiber-coupled measurements from this telescope were not useful during 
IHOP due to an improperly coated mirror. Because of this problem, most data products were processed to a 
minimum altitude of 300 meters; the lowest altitude to which the data from the large telescope were useful. 
The laser used was the same laser that was installed in the SRL prior to the TARFOX field campaign: a 
Continuum tripled Nd:YAG laser operating at 354.7 nm operating at 30 Hz and generating approximately 
9 W of output power. The fields of view of the large and small telescopes were approximately 0.3 and 1.2 
mrad, respectively. The narrow field-of-view of the large telescope coupled with narrow spectral filters in 
the water vapor and nitrogen channels (2.5 A and 2.8A, respectively) and the combined use of analog and 
photon counting data acquisition permitted the full intensity water vapor and nitrogen signals to be sensed 
throughout the diurnal cycle with no instrumental changes occurring during the measurements. The previ- 
ous configuration of both the SRL and CARL required attenuation filters during the daytime due to the high 
solar background. The addition of analog data acquisition electronics and narrow' spectral bandpass filters 
now permits the full intensity signal to be sampled diumally. This enabled the w'ater vapor mixing ratio to 
be quantified during the daytime at significantly higher temporal and spatial resolution than was previously 
possible. Table 1 shows the specifications of the SRL during IHOP and notes the components/measurements 
that were new for IHOP. 
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NASA/GSFC Scanning Raman Lidar Specifications for the IHOP field campaign 

Component 

Description 

Notes 


Laser 

Continuum Custom longcavity Nd:YAG 

converted to closely emulate model 


tripled. 30 Hz, 300 mj/pulse 

9030 

Telescope (high altitude) 

Starr Optical, Dall-Kirkham, 0.76 m, 
F/5.2 horizontally mounted 

variable FOV, but nomimally 0.3 milli- 
radians 

Telescope (low altitude) 

Custom Newtonian, 0.25 m, F2.5 

1 .2 milliradian field of view. New for 
IHOP. 


Spectrum Analyzer Custom designed using Barr dichroic 

beamsplitters and interference filters 


Measurement/ Filter - Large telescope 

Rotational Raman (352/1.2 newforlHOP 

nm) 

Rotational Raman (354/0.3 new for IHOP 

nm) 

Rayliegh-Mie (354.7/0.3 nm) 

Raman Nitrogen (386.7/0.3 

nm) 

Raman Liquid/Ice water new for IHOP 

(403/6.0 nm) 

Raman water vapor now tilt-tuned to center at 407.48 nm 

(407.6/0.25 nm) 


Measurement/ Filter - Small telescope 


Rayleigh-Mie (354.7/0.3 nm) 

fiber coupled 

new for IHOP 

Perpendicular polarization 

(354.7/0.3 nm) 

direct coupled 

new for IHOP 

Parallel polarization 

(354.7/0.3 nm) 

direct coupled 

new for IHOP 

Raman nitrogen (386.7/0.3 

nm) 

fiber coupled 

new for IHOP 

Raman liquid/ice water 

(403/6.0 nm) 

fiber coupled 

new for IHOP 

Raman water vapor 

(407.65/0.25 nm) 

now tilt-tuned to center at 407.48 nm 

new for IHOP 

Scanning Mirror 

Starr Optical, 1.1m x 0.8m flat, 50mm 
thickness 


Scan Motor Assembly 

Custom twin-axle, belt-drive design 
with 22-bit shaft encoder 


PMTs 

Hamamatsu R1924 (selected for high 
gain, high blue sensitivity and low dark 
counts) and R7400 


PMT Housings 

Products for Research 


Data Acquisition System 

Licel Analog + Photon Counting Tran- 
sient Recorders (7.5m spatial resolu- 
tion) 


Data Acquisition Computer 

PC running Windows 2000 



Table 1 : Hardware and optical configuration of the NASA/GSFC Scanning Raman Lidar during the IHOP 
field campaign 








5 Data Analysis Techniques 


Numerous atmospheric quantities were measured by the SRL during IHOP as shown in Table 1. For all of 
these measurements, the first step in the analysis is combining the AD and PC data into a single composite 

profile for each signal measured. That procedure will now be described. 

5.1 Combined Analog and Photon Counting Data - ’’Gluing” 

The Licel data acquisition electronics in use in the SRL measure a lidar signal simultaneously using 20 
MHz, 12-bit analog to digital converters and 250 MHz photon counters. Two separate data streams are 
maintained within the instruments that must be merged in some fashion to create a final output profile for 
each channel. The process of combining the analog and photon counting data together has come to be 
known as "gluing". After experimenting with various methods of calculating the conversion between an 
analog voltage and a photon countrate, the following procedure was used to determine mean conversion 
factors that are used to convert the analog data to a "virtual" countrate scale. 

The photon counting data are first corrected for photon pulse pile-up assuming the photon counting cir- 
cuitry behaves in a non-paralyzable fashion [Whiteman et. al., 1992] [Whiteman et. al., 2003a]. The back- 
ground is subtracted separately from the analog and photon counting records. Then, ordered pairs of analog 
and photon counting data are formed in a region of the signal where both are considered to be responding 
linearly and that avoids problem areas such as near the time when the laser fires. The low and high coun- 
trate thresholds that define this region are typically 1 and 10-20 MHz, respectively, for our Hama ma tsu 1924 
PMTs that possess a pulse width of 4-5 ns. The points in the profile that fall within these ranges and satisfy 
the other conditions are then formed as ordered pairs and a linear regression is performed on these points. 
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As an additional noise filter, the residuals between the actual data points and the best fit linear regression 
are then calculated and any points outside of 2 standard deviations from the mean are excluded. For the 
Rayleigh-Mie channels, additional noise filtering using 95% quantile regression was sometimes required. 
This regression technique is illustrated in the upper and middle plots of figure 1. In the upper plot, the 
regression of the analog and photon counting ordered pairs is shown with and without the application of the 
photon pulse pile-up correction. The slope (in MHz'mV) and offset (in MHz) of the regressions are given in 
the figure. Notice that the regression of the pulse pile-up corrected data exhibits a slope approximately 10% 
higher than the uncorrected data Also, the offset of the corrected data is much closer to zero. The middle 
plot of the figure is a blowup of the portion of the regression near the origin to permit the offsets to be 
seen more clearly. Since the background has already been subtracted separately from the analog and photon 
counting data, the value of the optimum resolving time can be deter min ed by studying the offset versus the 
resolving time used in the pile-up correction. The optimum resolving time is the one that yields an offset 
closest to zero. This technique can be contrasted with what was done previously to determine the op tim um 
resolving time using purely photon counting systems [Whiteman et. al., 1992] [Whiteman et al., 2003a] 
where a significantly attenuated version of the signal was used as the linear signal instead of the analog 
signal as used here. The technique illustrated here should be valid as long as the non-paralyzable equation 
properly describes the counting behavior of the electronics and the analog electronics are providing linear 
measurements. The non-paralyzable equation is certainly just an approximation of the counting behavior of 
the electronics so it is important to limit the maximum countrate of the data used to determine the optimum 
value of the resolving time. The analysis of the offset versus resolving time shown in the lower panel of 
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figure 1, performed using a countrate range of 1-20 MHz, shows that ~5 ns is the optimum value. For 
comparison, the RMS differences between the non-paralyzable model and the actual data points are also 
shown as a function of resolving time. There is a poorly defined minimum at approximately 2 ns. We do not 
consider this dete rmina tion of the resolving time to be reliable due to the very shallow minimum in the RMS 
curve. The minimum becomes better defined using this technique if higher countrates are permitted in the 
regression. However, considering that the non-paralyzable correction already amounts to ~1 1% at 20 MHz 
using 5 ns for the resolving time, determining the optimum resolving time using higher countrates will be 
more subject to errors introduced by the assumption that the electronics behave in a purely non-paralyzable 
fashion. Therefore, we choose to determine the optimum resolving time based on the zero crossing for the 
offset. 

5.1.1 Gluing coefficients and water vapor mixing ratio calibration 

Variations in the Raman lidar water vapor mixing ratio calibration can have important consequences for 
meteorological application of Raman lidar measurements. Therefore, it is important to study the behavior 
of the analog to photon counting conversion factors under differing conditions. In particular, the behavior of 
the slopes of the regressions as the background is increasing due to increased solar radiation is of particular 
interest. Data acquired on June 19-20, 2002 will now be used to illustrate the gluing process. 

The times series of fully calibrated water vapor mixing ratio is presented in figure 2a. As will be dis- 
cussed in section 5.2, the water vapor mixing ratio calculation is performed using the Raman water vapor 
and nitrogen measurements. Figures 2b and 2c therefore show which portions of the fully processed water 
vapor image shown in a) used AD (shown in white) and PC (black) data from the water vapor and nitro- 
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Figure 1 : Upper: two regressions of {analog, photon couting} data pairs are shown: before and after the 
correction for pulse pile up. Middle: the same data as in the upper plot but a blow up near the origin to 
show the intercept values more clearly. Lower: calculated offset of the linear regression versus resolving 
time value. The optimum resolving time of ~5 ns produces a near-zero offset. For comparison, the RMS of 
the residuals for each regression are also shown. This method of determining the optimum resolving time 
is less sensitive than the offset approach. 





gen signals. One can see that for daytime measurements only analog data were used whereas during the 
nighttime a combination of AD and PC are used. It is during the periods of reduced solar background when 
the gluing process can be used in the Raman channels. A time series of the gluing conversion factors is 
therefore presented in figure 2d for both the water vapor signal and the nitrogen signals where a transi- 
tion period between daytime and nighttime measurements has been chosen. At the beginning of the period 
shown, mean photon counting signals exceed 20 MHz and therefore a regression between AD and PC data 
was not performed. At approximately 0150 UTC (indicated as 2550 in the figure), the solar background 
decreased to the point where it was possible to perform the regression. But as the figure shows, the regres- 
sions that occur during the transition from day to night are characterized by low correlation coefficients (R 2 
multiplied by 10) and significantly changing slopes. Through experimentation, it was determined that good 
quality regressions resulted if the mean background countrate was below the minimum threshold frequency 
set to be 1 MHz in this example. The time at which this occurred (-2610 UTC) is indicated by the vertical 
line in the figure. The calculated slopes and correlation coefficients (R 2 ) tend to decrease quickly to the 
left of this line, which defines the region where the background points qualify for the regression. To the 
right of this line the mean H 2 O (indicated by m^HnO)) andN 2 (mfA^)) slopes are -8.4 ± 0.08 x 10 10 and 
9.5 ± 0.03 x 10 10 Hz/Volt (10 w Hz/Volt = 10 MHz/mV), respectively. The N 2 slope is quite constant 
beyond the transition point but, even though the variation in the H 2 O slope beyond the transition point is 
less than 1%, the slight tendency for the H 2 O slope to increase with time is still under investigation. 

Even though good quality regressions are obtained under most circumstances using the technique just 
described, the profile to profile conversion factor between AD and PC is subject to statistical fluctuations 
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Figure 2: First panel: water vapor mixing ratio variation during the period of June 19-20, 2002 in IHOP. 
Three wave regions are encircled. The major bore event was at -30.5 UT and 0.5 km. Second panel: 
diagram indicating the use of either analog (white)^ photon counting data (black) for the nitrogen data 
used in calculating the water vapor mixing shown in the top panel. Hurd panel: diagram indicating the use 
of either analog (white) or photon counting data (black) for the water vapor data used in calculating the 
water vapor mixing shown in the top panel. Fourth panel: time series of slopes of the regressions used to 
determine the coefficients to convert analog data to a virtual photon counting scale. 



that likely do not reflect variation in the gain of the electronics. Therefore, it is desirable to use a single 
conversion factor for each channel for an entire data record For the IHOP data processing, the mean values 
of the conversion factors obtained for a portion of nighttime data were thus used to glue the AD and PC 
data for those parts of the profile where the countrate exceeded 20 MHz. When the mean values of the PC 
background exceeded 1 MHz (e.g. prior to 2610 UTC in figure 2), the AD data (converted to countrate 

using the conversion factors) were used exclusively without gluing. 

5.2 Water vapor mixing ratio calculation 


The water vapor mixing ratio is defined as the mass of water vapor divided by the mass of dry air in 
a given volume. It is conserved in atmospheric processes that do not involve condensation or evapora- 


tion and thus serves as a tracer for air parcels. Also, the vertical profile of mixing ratio strongly influ- 
ences atmospheric stability. Water vapor mixing ratio can be calculated using the following expression 
[Whiteman et. al., 2003b] by taking the ratio of the Raman water vapor (407.5 nm) and nitrogen signals 
(386.7 nm) and accounting for the atmospheric differential transmission that occurs at the two different 
return wavelengths. 
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where w is the water vapor mixing ratio typically expressed in units of g kg -1 . Oh (r) /Oh (r) is the ratio 
of overlap functions for the N 2 (N) and H 2 O (H) channels, k is a proportionality constant determined by 
the ratio of molecular weights of water ( MWh 2 o ) and dry air {MW Dry Air) and the fraction of air comprised 
by molecular nitrogen (-0.78), Fh{T)/Fh{T) is the ratio of the temperature dependent functions for the 
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Raman jV 2 and H2O channels, where for example Fjj(T) is defined by 


Fh{T) = 




Fh{T) carries all the temperature dependence of the lidar equation for the water vapor channel. It contains 
the effects of any changes in the system transmission efficiency, £ (A) , for wavelengths other than A# 
within the passband A A//. £ (A//) is the transmission efficiency at A#. The notation do H {t) /dfl is used 
to indicate the total Raman backscatter cross section for water vapor at the stimulating wavelength. The 
product Fh{T) may be viewed as the effective molecular cross section that is consistent with the use 
of a monochromatic optical efficiency term, £ (Ajj) , in equation 1 . P (A*,r) = S (X x ,r) - B {X x ,r), 
where P ( Xx , r) is the background subtracted power in the Raman channel for species X, and 5 and B 
represent the signal and background, respectively. Thus, P (A jj, r) jP (Ajv, r) is the ratio of background- 
subtracted signals from the H2O and N2 channels, d °^ - is the ratio of the full Raman cross 

sections for No and HoO, £ (A^r) /£ (A#) is the ratio of N2 and H2O lidar channel efficiencies at their 
characteristic wavelengths of X x and A h and Ar (X H , A,v, r) is the differential atmospheric transmission 
that occurs at the Raman shifted wavelengths Xn and A h- The value of Ar is determined from: 


At (X h , A at, r) = exp J r) - a(X N , r')}dr'^ (4 

where r is the range to the volume of interest, and a is the extinction coefficient (units of inverse length). 


The equations used to quantify the random error in the water vapor mixing ratio, assuming Poisson 
statistics, are given below. 
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°%H + , °SN + °%N 

Rw ( Sh - Bh ) 2 (Sn ~ Bx) 2 

2 — (fta ~ Bh) 2 ( ct 5H ~P a BH , a |-N + a BN \ 

*- (Bn - Bn) 2 { (Sh - B H ) 2 + (S N - B N ) 2 ) 

In these equations, Rw represents the ratio of the background-subtracted water vapor and nitrogen lidar 
signals. In other words, Rw = (Sh — Bh) / (Sn — Bn) , where Sx and Bx are shorthand for S (Ax, r) 
and B (Ax , r ), respectively [Whiteman et. al., 2003b]. 

5.2.1 Overlap correction for the water vapor mixing ratio measurement 

Simple geometrical ray trace considerations indicate that in calculating a quantity from a ratio of two lidar 
channels that use a common field stop (such as the water vapor and nitrogen signals used in calculating the 
water vapor mixing ratio given by equation 1), the overlap functions On (r) and Oh ( r ) are equal and thus 
cancel. In real applications, however, this ratio does not equal unity in the near field and a residual overlap 
function must be determined and applied as a correction to the data. The existence of a residual overlap 
function appears to be due to the fact that as laser fight propagates from the near field to the far field of 
the telescope, fight received by the optical detection system spreads across optical components such as in- 
terference filters and photomultiplier tubes that may possess position-dependent efficiencies. These effects 
are greatest in the near field where modem radiosondes can provide high quality measurements of relative 
humidity, temperature and pressure. For the IHOP analysis, therefore, a mean residual overlap correction 
function was calculated by using 26 comparisons of SRL and Vaisala RS-80 radiosonde profiles of water 
vapor mixing ratio that occured during IHOP. The mean residual overlap correction function resulting from 
that calculation was then applied to all IHOP water vapor mixing ratio data as a part of the data reduction. 
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The residual overlap correction function was unity above an altitude of 750 meters and decreased to 0.94 
at an altitude of 300 meters, the minimum altitude of processed data for IHOP. Therefore the maximum 
correction produced by this residual overlap correction function was 6%. 

5.2.2 Temperature dependence of Raman scattering 

The temperature dependence of narrowband water vapor and nitrogen Raman measurements was also ac- 
counted for in the analysis of the IHOP water vapor mixing ratio data. An analysis of the temperature 
dependence of Raman water vapor and nitrogen scattering and its influence on the mixing ratio calculation 
indicates that 1) the effect is dominated by the temperature dependence of the water vapor spectrum and not 
the N 2 spectrum and 2) the net effect for likely filter configurations is to yield an apparent excess in water va- 
por concentrations that tends to increase with altitude [Whiteman et. al., 2003a] [Whiteman et. al., 2003b]. 
This latter effect is due primarily to a shift of intensity in the Raman spectrum of water vapor toward the 
band origin as temperatures decrease with altitude. 

To assess the magnitude of this effect for the IHOP measurements, laser output and interference filter 
transmission properties were carefully measured in order to apply the results of the earlier theoretical stud- 
ies. The Continuum laser frequency-doubled wavelength was measured using a Burleigh 4500 wavemeter 
and found to be 532.07 +/- 0.005 nm. Assuming errors are random and uncorrelated, this implies a fre- 
quency tripled wavelength of 354.71 +/- 0.003 nm. In addition, the water vapor and nitrogen interference 
filter bandpass characteristics were measured using a Thermo-Electron Nicolet 870 Fourier Transform Spec- 
trometer (FTS) operating at 0.5 cm' 1 resolution. The spectrometer was calibrated using a mercury lamp 
since the internal calibration based on a Helium Neon (HeNe) laser was found to not be reliable for measure- 
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ments in the near UV region of the spectrum (i.e. at shorter wavelengths than the HeNe calibration source). 
The combined error of the FTS measurement and mercury calibration did not exceed 0.001 nm. Figure 
3 shows the mercury lamp calibration corrected measurements of the water vapor interference filter trans- 
mission overlaid on the Raman scattering spectrum of water vapor calculated at 210K using the data from 
reference [Avila et. al., 1999]. The figure shows that the peak of the water vapor interference filter used for 
the IHOP water vapor measurements was shifted long of the peak in the Raman water vapor spectrum by 
approximately 0.05 nm. (Our comparison of traditional grating spectrometer measurements of filter cen- 
tral position and those obtained with the mercury lamp calibrated FTS instrument have revealed consistent 
differences in center line position of between 0.02 - 0.05 nm. The FTS measurements are considered more 
reliable due to the instrument’s inherent linear response, the repeatability of the mercury calibration of the 
FTS and, in a grating spectrometer, the difficulty of accounting for the sinusoidal variation of spectral posi- 
tion between calibration points.) The fact that the filter was positioned long of the peak of the water vapor 
spectrum introduced considerable temperature sensitivity to the measurement of water vapor over the range 
of temperatures present in the troposphere as will be shown later. It should be mentioned here that tilting 
the filter so as to operate at a shorter wavelength such as 407.45 would essentially eliminate the temperature 
dependence of water vapor with altitude. Since the time of the IHOP field campaign, this tilt-tuning has 
been implemented in the SRL. 

A study of the transmission properties of the SRL indicated that all significant variation in the lidar 
system efficiency over the passband of the Raman water vapor feature was confined to the water vapor 
interference filter itself. Therefore, the filter transmission measurements versus wavelength can be used 
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Wavelength (nm) 


Figure 3: Water vapor interference filter transmission measurements using a mercury-lamp calibrated 
Fourier Transform Spectrometer operating at 0.5 cm -1 resolution. For reference, the Raman water va- 
por spectrum simulated at 270K is also plotted. The shift in the filter from the peak of the Raman spectrum 
implies significant temperature sensitivity for the water vapor measurement. 
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to indicate the variation of the lidar system efficiency over the water vapor passband. Similar analysis of 
the Raman nitrogen measurements were made and combined to generate the mean temperature correction 
vector, Fff [ T (r)] /Fjj [T (r)] from equation 1, that was applied to all the SRL water vapor mixing ratio 
profiles during the reduction of the IHOP data. This correction vector is shown in figure 4. The mean 
IHOP temperature correction vector varies by approximately 10% from the surface to an altitude of 14 km. 
This implies that, for the configuration of the SRL during IHOP, the correction to the water vapor mixing 
ratio due the temperature dependence of Raman scattering is approximately 10% in the troposphere. For 
comparison, the temperature correction vector calculated for the same filter characteristics but assuming 
the U.S. Standard Atmosphere temperature profile is also shown. As this comparison implies, the residual 
error due to using a single temperature correction vector for all of the IHOP data is estimated to be 1% or 
less. The temperature correction vector for the same set of water vapor and nitrogen filters, but if the water 
vapor filter had been tilt-tuned to 407.45 nm, is also shown in the figure. Notice that for a filter position of 
407.45 nm, the values of the ratio Fjv \T (r)] /Fjj [ T (r)] are lower for this nearly temperature insensitive 
configuration indicating that a larger fraction of the Raman water vapor cross section is being transmitted. 
Figure 4 demonstrates that for the configuration of IHOP a significant temperature- (and thus altitude-) 
dependent correction was required for the water vapor mixing ratio calculation. However, the figure also 
shows that, even for a narrow water vapor interference filter such as used in the SRL for IHOP (2.5 A), the 
temperature sensitivity can essentially be eliminated with proper tuning of system parameters if accurate 
measurements of laser output wavelength and lidar system transmission characteristics are available. 
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Figure 4: Temperature correction vectors Fjy [T (r)] /Fjj [T (r)] for the water vapor mixing ratio calcula- 
tion. The vector applied during the IHOP analysis is shown along with the vector corresonding to the U.S. 
Standard Atmosphere assuming the same experimental configuration as in IHOP. The curve laveled 407.45 
is the temperature vector that would result if the water vapor filter had been positioned at 407.45 nm through 
tilt tuning. This latter configuration is nearly temperature insensitive throughout the troposphere. 
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5.2.3 Water vapor mixing ratio calibration 


A first principles Raman water vapor lidar calibration can be performed with standard optical labora- 
tory procedures. However, the ratio of the Raman cross sections for water vapor and nitrogen is not 
known to better than 10%, which implies that the absolute accuracy of such an effort will exceed 10% 
[Vaughan et. al., 1988] [Sherlock et. al., 1999]. Calibration by comparison with other water vapor sensors, 
such as research grade radiosonde, microwave radiometer or GPS, is thus the standard within the Raman 
water vapor lidar community and the general technique that was used to calibrate the SRL water vapor 
measurements during IHOP. 

In particular, the SRL water vapor mixing ratio measurements during IHOP were calibrated by com- 
paring the integral of the lidar mixing ratio profile with the total precipitable water (PW) derived from a 
SuomiNet GPS [Ware et. al., 2000] system mounted on the SRL. During the daytime, the SRL water vapor 
profile extended usefully only to the top of the boundary layer, thus leaving a significant fraction of the 
total precipitable water unmeasured. Therefore to calculate PW from the lidar during the daytime, SRL 
calibrations with respect to GPS were limited to radiosonde launch times and the radiosonde profile was 
normalized to the lidar and used to extend the lidar mixing ratio profile upward beyond the point at which 
25% random error was present in the lidar data. To account for the portion of the mixing ratio profile 
near the ground not measured by the lidar, a linear interpolation was performed between the lowest altitude 
measured by the lidar (typically 300m) to the ground point measured by a Paroscientific Met3A station as- 
sociated with the SuomiNet GPS system. This composite profile was then integrated to yield the SRL total 
precipitable water for daytime comparisons with GPS. For nighttime comparisons, it was not necessary to 
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extend the profile upward using radiosonde since the lidar water vapor data were of sufficient quality to 
permit —99% of the PW to be quantified based on standard atmosphere concentrations of water vapor above 
the maximum lidar height measured. However, the downward extension to the ground was still necessary. 

Except for three dates during IHOP (June 14, 1 7, 1 8), the water vapor data were processed using a single, 
height independent calibration constant determined from the GPS calibration procedure just described. On 
these days, the lidar calibration differed by 10-15% from the mean calibration value perhaps due to acci- 
dental changes in system operating parameters such as PMT high voltage. For these cases the calibration 
constant used was the one determined on that day instead of the mean value. The standard deviation of the 
calibration constant for the IHOP water vapor mixing ratio data, including the effects of these anomalous 
days, was approximately 6%. When using only daytime measurements it was approximately 6.5% and 4.5% 
when using only nighttime measurements. The smaller standard deviation of the calibration constant for 
nighttime measurements is thought to be due to greater atmospheric horizontal homogeneity at night and 
thus the better agreement between the profile of water vapor measured over the lidar site and the volume 
average measurement of the GPS. The daytime and nighttime calibration constants agreed to within -1% 
implying that there was no significant difference in the lidar water vapor calibration constant due to diurnal 
effects. 

GPS as the calibration source for Raman water vapor lidar The U. S. DOE Atmospheric Radiation 
Measurements (ARM) program calibrates its Raman lidar, CARL, in a similar fashion as just described for 
the SRL except that the CARL PW is compared to the total column water vapor measured by microwave 
radiometer [Turner et. al., 2002]. Research done within the DOE ARM program indicates that carefully 
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calibrated and analyzed microwave radiometer (MWR) data possess an absolute accuracy of approximately 
3-4%. This makes it an excellent calibration standard for atmospheric research. A deployment of the SRL 
to the DOE Southern Great Plains site in Oklahoma in 2003 for the AIRS Water Vapor Experiment-Ground 
(AWEX-G) [?] permitted a careful comparison of the ARM microwave radiometer and the same SuomiNet 
GPS system that accompanied the SRL into the field for LHOP.The GPS measures over a much larger 
volume than the MWR and therefore individual comparisons can show considerable disagreement under 
conditions of spatial non-homogeneity in the atmosphere. Line of site comparisons of the two instruments 
have been performed to address these differences and have shown excellent agreement [Braun et. al., 2003]. 
During AWEX-G, we performed an extended comparison of 30-minute average GPS and MWR vertical 
precipitable water measurements in order to minimize the effects of short term spatial inhomogeneities. 
The results of that comparison are shown in figure 5. 

The plot on the left of the figure shows that the slope of the regression between GPS and MWR was 
essentially unity but, due to an offset of ~0.4 mm, the PW from GPS was on average 2.4% higher than 
MWR This overall agreement of the two sensors is within the uncertainty of the MWR supporting the use 
of the GPS system as an independent source for calibration. The plot on the right of figure 5 presents the 
comparison of MWR and GPS PW hour by hour. The overall moist bias of the GPS with respect to MWR 
that was indicated by the regression can be seen throughout the diurnal period. However, the GPS shows 
a slight tendency to become moister during the daytime period from —1200 - 2400 UTC. The source of 
this small diurnal bias is under investigation but it is not thought to be due to variations in the MWR since 
comparisons of AERI and MWR scaled radiosondes have indicated no significant diurnal bias in the MWR 
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SuomiNetSA35 vs SGP MWR 




Figure 5: Left: An analysis of the precipitable water measurements from the SuomiNet GPS that accompa- 
nied the SRL to IHOP and served as the source of the water vapor calibration and the U.S. DOE microwave 
radiometer at their Southern Great Plains research site. There is a mean offset of -2.4% which is within 
the absolute uncertainty of the MWR. The best fit regression is shown with a dashed line. The solid line 
indicates perfect 1-1 correlation. Right: The same data as plotted on the left except this time divided by 
time of day. A small diurnal bias between GPS and MWR is observed 
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[Turner et. al., 2002], This analysis indicates that the GPS calibration agrees well with the DOE MWR 
and offers similar absolute accuracy as a source for Raman water vapor lidar calibration provided sufficient 

statistics are accumulated to reduce the effects of spatial non-homogeneities in the atmosphere. 

5.3 Aerosol scattering ratio and backscatter coefficient 


The aerosol scattering ratio, 71, is de fin ed as the ratio of the total (molecular and particle) backscatter 
coefficient divided by the molecular backscatter coefficient. It can be calculated from the ratio of the 
received power in the Rayleigh-Mie and Raman vibrational N 2 channels. The equations for calculating the 
aerosol scattering ratio, including the effects of the temperature dependence of rotational and vibrational 
scattering are [Whiteman et. al., 2003b] 


n(\ L ,r)-l = C' N (X L , r) F N (T (r)) f At (X n , X l , r) - F R (T (r)) (7) 


C* N {X L ,r) = C N (X L ) 


P(AX n ,t) 
On (r)£(XN) 

Or (r) £ (Xi) 


( 8 ) 


where P (AA r, r) and P (AA/y, r) are the background-subtracted received power in the channels measur- 
ing the Rayleigh-Mie signal in the spectral band AA r and the vibrational N 2 Raman-shifted in the spectral 
band AA^r- The differential transmission. At (Ajv, Ax, r) , is the ratio of atmospheric transmission at the 


two wavelengths, Ajv and A x and is calculated using an equation similar to 4. The lidar channel optical 
efficiencies are expressed as f (A^r) and £ (Ax). The lidar system overlap functions are given by On ( r ) 
and Or ( r ) . The calibration constant Cjv has been evaluated at the laser wavelength of ~355 nm using 
the values of the Raman vibrational and Rayleigh differential scattering cross-sections, dc (it) /dQ, and 
dajnoi ( tt ) /dQ, respectively [Measures 1984], The effect of the temperature dependence of the Raman 
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scattering on the Rayleigh-Mie and Raman nitrogen signals is contained in the two terms, Fr(T (r)) and 
Fk(T (r)), respectively, which are calculated using an equation similar to 3. The aerosol backscatter coef- 
ficient, /? Qer (Ax, r), can be evaluated from the aerosol scattering ratio as follows 

/^(A x, r) = (Al, r) (K (Ax, r) - 1) 

where /?™ oi (Ax, z) is the Rayleigh backscatter coefficient at the laser wavelength calculated using density 
measurements from a radiosonde. 

5.4 Cirrus cloud optical depth and layer mean extinction to backscatter ratio 

Cirrus cloud optical depth can be calculated from a Raman lidar measurement of molecular nitrogen, which, 
if properly performed, shows only attenuation due to the presence of the cloud. The amount of this atten- 
uation can be converted to optical depth once the atmospheric density is known. The single scattering 
equation which yields two-way optical depth is obtained by integrating the equation for aerosol extinction 
[Whiteman et. al., 2003a] and can be written as: 


[ [a(Ax, r) + o(Aas r)j dr 
Jrj 

, /0^(r- 2 )i r i v(T(r 2 ))iV J v(r2)r2p(A AS r 1 )^ f r \ ,, N 

= ln ( o,(n)Wr(n))iV.(r,) r |P(l„.r 2 ) j-X, + * 

in ( - f MltO + MUr)] dr 

f, r 2) ) Jr , 


( 10 ) 


V Nn ( r l) r 2 P (A N, 

where r\ is below the cloud, r 2 is above the cloud. Ax is the laser wavelength (354.7 nm), A_v is the 
wavelength of the Raman nitrogen signal (386.7 nm), a(A x , r) is the cloud extinction coefficient as a 
function of wavelength and range, On (r) is the overlap function for the nitrogen channel evaluated at 
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range r, F^(T (r)) is the temperature dependent factor for the Raman nitrogen measurement, Njy (r) is 
the number density of atmospheric nitrogen (using the full atmospheric number density yields equivalent 
results) as a function of range, P(Ajv,r) is the background-subtracted Raman lidar nitrogen signal and 
c*moi(A x , T ') is the extinction coefficient due to molecular scattering obtained from radiosonde data. At 
typical cirrus altitudes, the ratios of the overlap and temperature dependent factors are nearly unity and 
can be ignored as shown in the final form of equation 10, which is the same result that would be obtained 
by integrating the traditional aerosol extinction equation [Ansmann et. al., 1992]. Assuming no multiple 
scattering and that cirrus cloud extinction is wavelength insensitive between A l and Ajv, the optical depth 
of the cloud is simply one-half of equation 10. 


A modified approach to the evaluation of equation 10 was used here which implements an iterative 
procedure that corrects for the influence of multiple scattering [Whiteman et. al., 2001a] using a Gaussian 
approximation technique [Eloranta, 1998] and calculates a layer mean particle radius in the process. The 
integrated backscatter is determined by integrating the profile of cirrus cloud backscatter coefficient, which 
has been shown to be essentially insensitive to multiple scattering [Wandinger, 1998], The layer mean 
extinction to backscatter ratio, also known as the lidar ratio, is then just the ratio of the cloud optical depth 

and the integrated backscatter coefficient. 

5.5 Aerosol depolarization 


The ability to calculate scattering ratios using the Raman lidar technique permits both the volume and 
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particle depolarization ratios 6 vo i and S paT to be calculated [Behrendt and Nakamura, 2002] as follows. 


Svoi (z) 
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where j3 refers to the backscatter coefficients for either molecules (mol) or particles (par ) in either the 
perpendicular (X) or parallel (||) directions. Rx is the scattering ratio for the perpendicular polarization 
signal and is the scattering ratio for the parallel polarized signal. The equation for the scattering ratio 
used here is given in 7 and fully accounts for both the temperature dependence of Raman scattering as 
described in section 5.3. 

The aerosol depolarization measurements were calibrated by first determining the relative gain of the 
parallel and perpendicular channels by observing a fully depolarized source, in this case thick overcast skies. 
Using this value of the relative gain in the channels, the volume depolarization ratio in clear air was found 
to be -1.7% as opposed to the theoretical value of molecular depolarization of -0.5% that corresponds 
to our -0.5A bandpass filters. The additional depolarization was attributed to misalignment between the 
planes of polarization of the outgoing laser and the polarization analyzer [Reichardt et. al., 2003], In order 
to compensate for the resulting cross-talk between the parallel and perpendicular channels, an approximate 
correction technique was used here. The volume depolarization ratio was calculated assuming variable 
amounts of cross talk between the parallel and perpendicular channels. The amount of assumed cross talk 
that yielded the theoretical value of 0.5% depolarization in clear air was used to analyze the depolarization 
data presented in part II. This amount of cross talk corresponded to an angular difference of —3 degrees 
between the laser polarization and the position of the polarization analyzer. Another correction that was 
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applied to the depolarization data was to account for the variation in the molecular depolarization as a 
function of temperature [Behrendt and Nakamura, 2002]. The value of molecular depolarization for the 
filters in use and the range of temperatures observed during IHOP was found to vary as a function of 
temperature from approximately 0.46 - 0.52%. 

6 Summary 

The NASA/GSFC participated in the first International H 2 O Project in May-June 2002. A new configu- 
ration of the SRL enabled measurements of water vapor, aerosol backscatter, extinction, depolarization, 
liquid water, ice water and rotational Raman temperature. In particular the combination of narrow field of 
view, narrow spectral bandpass, relatively large pulse energy Nd:YAG laser and both analog and photon 
counting detection electronics permitted a single configuration to be used for both daytime and nighttime 
measurements of water vapor mixing ratio. The improved water vapor measurement capability permits 
convective boundary layer processes to be studied throughout the diurnal cycle. In this first of two parts, the 
method of joining the analog and photon counting data, referred to as "gluing”, was described for the case 
of water vapor measurements during a period that included a large number of atmospheric waves. For the 
first time, the temperature dependence of Raman scattering was accounted for in the calibration of a Raman 
water vapor and aerosol lidar and shown, for the configuration of the SRL during IHOP, to effect the upper 
tropospheric water vapor measurements by approximately 10%. A study of the calibration of water vapor 
mixing ratio indicated a stable calibration constant to within 1% between daytime and nighttime measure- 
ments. The measurements presented here were all provided using a single output wavelength of 354.7 nm 
demonstrating the feasibility of offering all of these measurements in an automated, eye-safe Raman lidar 


37 



system. In part two [Whiteman et. al., 2005b], comparisons of SRL measurements with other sensors are 

presented along with daytime and nighttime case studies. 
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